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Abstract 

A new method is developed which permits the reconstruction of the surface- 
density distribution in the galactic disk of finite radius from an arbitrary 
smooth distribution of the angular velocity via two simple quadratures. The 
existence of upper limits for disk’s mass and radius during the analytic con¬ 
tinuation of rotation curves into the hidden (non-radiating) part of the disk 
is demonstrated. 
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The rotation curves of spiral galaxies supply astronomers with the direct information 
about distribution of the gravitational potential in galactic disks. The measurements of 
Hq, hi and CO lines allowed to construct rotation curves for our and nearby galaxies (for 
early reviews cf. Refs. [1-3]). The advanced modern techniques of observation and data 
processing permitted to extend rotation curves at the distances far beyond the galactic 
center and discover the flat rotation curves [3-7]. The hrst galactic disk models were based 
on the concept of oblate shells (homoeoids) between confocal spheroids [8-10]. Brandt 
0, Brandt and Belton |l^ tended the eccentricity of confocal spheroids to zero, thus 
obtaining the disks covering the whole plane ^ = 0. They proposed a special rotation 
curve possessing Keplerian asymptotics for M31 and NGC5055, and found numerically the 
corresponding surface-density distribution. Another approach was developed by Toomre |T3 


who considered a new class of asymptotically Keplerian rotation curves and found, for inhnite 
disks, the corresponding surface-density distributions. The problem of reconstruction of the 
surface density from the rotation curve in the case of inhnite disks containing a black hole 
at the center was solved in our work []^. Mestel |K] pointed out that the highly battened 
galaxies and nebulas should be modeled by self-gravitating disks of finite radius R. He 
discussed two limiting cases: a rigidly rotating disk and a model with a constant velocity. 
In a different context some exact models of self-gravitating disks with outer boundary were 
found by Morgan and Morgan [^ . 


In the present paper we propose a new method to the problem of reconstruction of the 
surface-density distribution in a self-gravitating disk of hnite extension from an arbitrary 
smooth distribution of the angular velocity, and it permits to obtain the general solution 
of this problem in terms of two simple quadratures. We shall also demonstrate that if a 
rotation curve is analytically continued into the hidden (non-radiating) part of the disk, 
then the disk’s mass can achieve its maximal value at a finite radius, and not necessarily at 
inhnity. Like it is done in the classical papers [10-13], we shall assume an idealized model 
for the highly battened galactic disks in which the inbuence of viscosity, pressure, non-disk 
component and all sort of non-circular motion of matter will be neglected. We will not 
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separate the gravitational fields of the dark (non-radiating) and radiating matter, firstly, 
because they both exert the same influence on rotation curves, and, secondly, because no 
reliable mechanism for such separation is yet known [17, 4, 18]. 

Our approach to the Newtonian self-gravitating disk problem is based on a new repre¬ 
sentation for the gravitational potential of a disk with the boundary radius R which we take 
in the form 

= \Jp^ + z^ ( 1 ) 

where the auxiliary space variables p, z are related to the genuine cylindrical coordinates 
p, 2 ; via the formulas 

p = R^p/r‘^, z = R^z/r‘^, = p^ + z^ , ( 2 ) 


and the function $ is defined in the following way: 

TT 

f de 

0 R 

q;(s) being a real function which has the meaning of the density of sources distribution. 
Due to the specific analytical properties of the logarithmic function which turn out to be 


s — pcosOy + z^ 



most advantageous for fulfilling the Laplace equation in the region exterior to disk’s sources, 
hereafter we shall be able to relate the function q;(s) to the surface density (t(p) of the disk 
and to the distribution of the angular velocity a;(p) via simple integral formulas. 

The potential ^{p,z) which, as can be easily seen from (1), is connected with $ by 
Thomson’s transform, satisfies Laplace’s equation everywhere except on the disk. On the 
circle z = d, p < R the normal derivative d^/dz experiences a jump. To show this, let us 
tend to zero. Then, by virtue of (2), we have 


(9<h p^ ,. <9$ 

hm = —- hm , 

z^+o dz R? i^+o dz 


while from (3) follows that 


lim ^ = lim ^ 

2^+0 dz z^+o 271^ 


-z —^-7 -rz H- ^ „ -zz I a(s) ds . 

s — p cos U — iz s — p cos U + iz J 
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Due to the Sokhotsky-Plemelj formula for the Cauchy-type integrals, we obtain from (5) 


0 < p < i? 


\ I f a(s)ds . 
li^+o - / , , p> R, 

TT J VP — 

R ^ ^ 


or, returning to the coordinates p, 2 :: 


(9$ —^ 

—— = —27rG(T(p) = { 

0 , 


Q;(s)ds 
— p^s 


= , 0 < p < i? 


p> R, 


where (j{p) is the surface density of the disk. 

The condition of balance of the gravitational and centrifugal forces for zero-pressure 
matter of the disk at 2 ; = 0, p < i?, i.e.. 


9$ 

z^O 


-\- uj p — 0 , 


can be represented in the form 




whence it follows that 


From (3) we obtain the form of the left-hand side of (10) in the limit z ^ 0: 

dp ^ [I I p's‘-p‘) 

In the subsequent formulas it is convenient to pass to the non-dimensional variables 


t = R^/s^ , a; = i?Vp^ = pV^^ 


In the new variables, formula (7) for the surface-density distribution assumes the form 


aix = - 


R f a{t) dt 
Att'^Gx J t\/t — X 
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(note that a(x) is not singular at a: = 0 because of the condition (15) below), so that we are 
able now to calculate the total mass of the disk: 

H 1 

rlT r mlt \ (it 

(14) 


M = 2nJa(p)pdp=-^I^J 


dx 

X 


a{t) dt 


0 0 
For M to be a hnite quantity, it is necessary that 


ty/t — X 


f a(t) dt f a(t) df 


= 0 . 


(15) 


0 


From (14) and (15) a simpler formula for the mass of the disk is readily obtainable: 

1 


M = 


f a(t) 


4:71G J RR 

0 


In t dt. 


( 16 ) 


Taking into account (11) and (15), the balance equation (10) yields the integral equation 
of the Abel-type with respect to a(t): 


a{t) dt 


27r J RR\/x — t 


X 


/ dt 

Vt 


+ Ao, 


(17) 


0 ' 0 

where the integration constant Aq should be found from the condition (15). 

The solution of equation (17) can be written in the form (cf. Ref. |^, p. 318) 


a[x) 

2x^R 




An 


dx J \fx 
0 


X 


dt 

y/iy/x — t 


(18) 

0 ■ -0 

In the general case the surface density has the following asymptotic behavior in the 

neighborhood of the outer rim x —»■ 1: 

2R f {ujH^^'^)[y/idt 


C 

cr(x) Ay/l - X, 


y/T^ 


(19) 


Formulas (13), (18) provide the general solution to the problem under consideration in 
terms of two successive quadratures: for a given distribution uj{p), one should first find the 
function a{p) with the aid of the formula (18) and then calculate the function a{p) from the 
integral (13). The total mass of the disk is dehned by (16). 

Let us illustrate the use of the general formulas with some concrete examples. Among 
numerous classes of smooth rotation curves which can be treated analytically with the aid 
of our formulas we choose the physically meaningful velocity distributions of the type 
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( 20 ) 


o;^ = j (6^ + ojqi b = const. 


This particular curve was proposed by Brandt |]^ in the case of infinite disks and it 
possesses the Keplerian asymptotics at p —> cxo; its parameter b multiplied by a/2 has the 
meaning of the radius of the circumference at which the linear velocity has its maximal value, 
while ujq is the angular velocity at the center of the disk. A wide family of asymptotically 
Keplerian rotation curves having the form 

n—1 .. 

( 21 ) 


= 




d 


n\ 


dh^ 


6(62 + ^ 2 ) 3/2 ’ 


n = 1,2, 


and containing (20) as the particular n = 1 specialization was studied by Toomre [0; the 
corresponding surface-density distributions (in the case of infinite disks) can be shown to 
be given by the formula 

62^+2 (2n- 1)!! 


(^nip) = 


( 22 ) 


2^n\ 7rG(62 + p2y2n+l)/2 

In what follows we shall obtain the analog of the expression (22) in the case of finite 
disks. For this purpose it is advantageous to rewrite the rotation curves (21) in terms of the 
dimensionless variable x = and parameter a = 6^/i?^: 

1 


= AO 

n 


,2 ^n+1 
u/n CL 


n\ 


da 


n—1 


(23) 


y/a{x + a)^/2 

(all these a;„ take the same value at a; = 0: a;„(0) = cuo). 

Let ns first consider the case n = 1. Setting n = 1 in (23) and substituting the resulting 
ujf into (18), we obtain 


afix) 

2x^R 


— UJc 


(y/x{3a + 


xja 


2 a 


\ (a + xfi 1 + 


(24) 


The corresponding formula for the surface-density distribution is then obtainable from 
(13), and it has the form 


aPx = 


a;2 aR%S 


^(1 — x)(a + x) + (1 + a)ArctanY/ (1 — x)/(a + x) 
^ = (l + a)(a + (r)3/2 ' 
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Due to the linearity of the problem, the surface-density distributions corresponding to 


the rotation curves (23) will be dehned in the general case by the expression 

,2^n+i^ / a \ n-l 


an{x) = 


lUq a 


nl 


da 


S. 


(26) 


In the limit R —> oo, formula (26) reduces to Toomre’s expression (22). In Fig. 1 we 
have plotted the surface-density distribution (25) as the function of p/h for different values 
of the ratio iP'/R?. 

The total mass of the disk dehned by the surface density (26) is given by the formula 

n—1 / 


,, 2ujI 


da 


11 1 ^ 

-1—^Arctan^^ 

1 + a \/a Ja, 


(27) 


n\ tiG 

and it can be readily found from (16). 

It is worthwhile mentioning that the rotation curves can be also approximated by linear 
combinations of different functions o;^ dehned by (21): 

N N 


00 




(28) 


n=l n=l 

and they will all possess the limit of the rigid-body rotation with a given angular velocity ujq 
near the center and the Keplerian asymptotics at large p. The corresponding surface-density 
distributions and total masses will be given by the formulas 

N N 


^ H c„an(x ), M ^ y CnMn{x ), 


(29) 


n=l n=l 

and in the neighborhood of the outer rim the surface density will have the following asymp¬ 
totics: 


a{x) 


oOqR^/1 — X 

^2^ 


N 


S '" ^“+1 


(30) 


An interesting important question which we would like to discuss in relation with our 
formalism is the possible estimation of upper limits of the galactic masses and radii of 
outer rims since, as is well known, real galactic disks can extend at much larger distances 


than might be thought from the observations of hot gas [10|. Here we shall restrict our 
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consideration to only the case of analytic continnation of the rotation cnrves into the region 
p > Ry where Ry is a known radins of the visible part of the galactic disk. 

Let an arbitrary smooth distribntion a;^(p) be known from observational data for 0 < 
p < Ry. Then the corresponding a{p) shonld be constrncted via the formnlas (13) and 
(18) at the snpposition that the real radins R of the disk is an nnknown parameter to be 
determined. The npper limit Rmax for the radius of the disk can be found from the condition 
(dcr/dp)p = 0 or A = 0 (see formula (19)), and the value of Rmax thus obtained should 
be subsequently substituted into the formula (16) to get the maximal value of the mass. A 
different, though equivalent, way of obtaining Rmax is to hnd it from the condition of mass’ 
maximum. Hence, for the real radius R and real mass M of the galactic disk we get the 
inequalities Ry < R < Rmax and M < Mmax- Mention, that for some rotation curves Rmax 
may have finite values, but it can also happen that Rmax = oo (like in the case of the pure 
Toomre’s angular-velocity distributions o;^, so that only a linear superposition of different 
o;^ leads to hnite Rmax)- 

The existence of maximal values for disks’ masses at hnite radii for the analytically 
continued rotation curves constitutes a new physical phenomenon (very similar to the static 
criteria of stability for the barotropic models of stars [^]) which is absent in inhnite disks. 
For instance, in the case of the angular velocity distribution (28) the condition A = 0 yields 
the algebraic equation 



As can be easily seen, for all > 0, n = 1, A^, this equation has the unique non-negative 
root a = 0 corresponding to R = oo. But already for ci < 0, q > 0, i = 2, N, it will have an 
additional, positive root a = a* > 0. In the interval 0 < a < a* the surface density turns out 
to be negative in the vicinity of the outer rim x = 1, so expressions (29) are only physically 
meaningful for a > a*, the boundary value of the disk radius being Rmax = h/y/cR. 

Let us illustrate the above said with the simplest linear combination of the curves (20) 
possessing the Keplerian asymptotics: 
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(32) 


o 1-2 ^/-i \2 2 

u = —-—+ -(1 — q)iLi2 = uJq 

0<g<l, 0<a<0.5. 


\a + xj \a + xj 


\ 5 / 2 ' 


The corresponding surface-density distribution is given by the formula 

<^( x ) = '*'^3 b i ( i ) + |(1 - 9 ) ff 2 ( i ), 


(33) 


where ai and a 2 are dehned by (25). 

The derivative (da/dp)p=R will be equal to zero when —1 + 3a + 4g = 0, whence one 
obtains that a is equal to (1 — 4g)/3. 

On the other hand, the mass of the disk with the angular-velocity distribution (32) is 
equal to 


M = 


2uJnb' 


2u3 


TT 


y/a{2 — 5q — 3aq) 


+ gArctan 


(34) 


3(l + a)2 

and, as function of a for hxed this expression achieves its maximum at a* = (1 — 4g)/3. 

Therefore, we have arrived at the following important conclusion: 

When 0 < g < 0.25, the radius of the disk cannot be greater than b^3/{l — 4g). For 
0.25 < g < 1, as can be seen from (34), the mass of the disk becomes a monotonically 
decreasing function of a, so that its maximum is achieved in the limit R oo. The 
dependence of Mmax on g is shown in Fig. 2. This example reveals in particular the difficulties 
which one may encounter while approximating a real rotation curve by superposing standard 
rotation curves in the case of inhnite disks; these difficulties are absent in the case of disks 
of hnite extension. 

This work has been supported by Projects 34222-E and 38495-E from CONACyT of 
Mexico. N.R.S. acknowledges hnancial support from SNI-CONACyT, Mexico. 
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FIGURES 
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FIG. 1. The form of the surface density distribution (25) for several particular values of a 
(a = 0,0.16,0.25,0.35,0.5). 
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FIG. 2. The dependence of the maximally possible mass of the disk on the parameter q in the 
example given by formulas (32)-(34). 
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